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ABSTRACT 

We report the search for low-mass X-ray binaries in quiescence (qLMXBs) in the 
globular cluster (GC) NGC 6304 using XMM-Newton observations. We present the 
spectral analysis leading to the identification of three candidate qLMXBs within the 
field of this GC, each consistent with the X-ray spectral properties of previously iden- 
tified qLMXBs in the field and in other GCs - specifically, with a hydrogen atmo- 
sphere neutron star with radius between 5-20 km. One (XMMU 171433—292747, 
with .Roo = 11.61^6 (D/5.97kpc)km and kT e{i = 122+fjj; eV is located within one 
core radius (r c ) of the centre of NGC 6304. This candidate also presents a spectral 
power-law component contributing 49 per cent of the 0.5-10 keV flux. A second one 
(XMMU 171411-293159 with = 10.7±| f (D/5.97 kpc) km and kT cS = 115+ 2 geV) 
is found well outside the optical core (at ~ 32 r c ) but still within the tidal radius. From 
spatial coincidence, we identify a bright 2MASS infrared counterpart which, at the dis- 
tance of NGC 6304, seems to be a post-asymptotic giant branch star. The third qLMXB 
(XMMU 171421-292917 with = 23+^ (D/5.97 kpc) km and kT eS = 70+^ eV) is 
a low signal-to-noise candidate for which we also identify from spatial coincidence a 
bright 2MASS infrared counterpart, with 99.916 per cent confidence. Three qLMXBs 
from this GC is marginally consistent with that expected from the encounter rate of 
NGC 6304. We also report a low signal-to-noise source with an unusually hard photon 
index (a = — 2.0j^ 2 'i)- Finally, we present an updated catalogue of the X-ray sources 
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lying in the field of NGC 6304, and compare this with the previous catalogue compiled 
from ROSAT observations. 

Subject headings: Stars: neutron, X-ray: binaries, Globular clusters: individual: NGC 6304 



1. Introduction 



Transiently accreting low-mass X-ray binaries in quiescence (qLMXBs) produce useful con- 
straints on physical models of the interiors of neutron stars. They were first discovered during 
the post-outburst periods of transient LMXBs Cen X-4 and Aql X-l, in which their faint emis- 
sion (Lx ~ 10 32 — 10 33 erg s _1 ), interpreted as thermal black-body emission, had emission areas 
a factor of 100 smaller than implied by the projected area of a 1 km neutron star (NS). Further 
observations of C en X-4 and new observations of 4U 2129+47 (jGarcial Il994l ) and 4U 1608-522 
(jAsai et al.lll996bl ) confirmed these observational properties of the class of qLMXBs. At that time, 
it was suggested that some low- level mass accretion on to the compact object was required to 
power the observed luminosity, although wha t set the low-mass accretion rate (M) was unclear 
dVan Paradiis et alJfl987l : IVerbunt et aD 1 19941 ). 



An alternate explanation to low-M accretion has become the dominant explanation for the 
quiescent emission of qLMXBs (lBrown et al.lll998l . BBR98 hereafter). An up to date summary 
of the theoretical developments of deep crustal heating and the related observations is provided 
here. In this different interpretation of the quiescent thermal emission, the luminosity is provided 
not by accretion on to the compact object, but by heat deposited in the neutron star crust by 



pressure-sensitive nuclear reac tions. These reactions (jHaensel fc Zduniklll990l . 



2007 



2003 



Gupta et al 



Haensel h Zdunikl 120081 . HZ08 hereafter), referred to as deep crustal heating (DCH), occur 



as matter is piled at the top of the neutron star surface, forcing the column below to greater 
density and electron Fermi energies. For example, beginning with pure 56 Fe at the top of the crust, 
the increasing electron Fermi energy induces electron capture at a density of 1.494 xlO 9 gem -3 ; 
because of the odd-even staggering of nuclear masses, a second electron capture follows to produce 
56 Cr and releases 0.041 MeV per accreted nucleon (HZ08). At a density 1.114xl0 10 gem -3 the 
electron Fermi energy is again above threshold for capture on to 56 Cr, and 56 Ti is produced with an 
energy deposition of 0.036 MeV per accreted nucleon. This process continues, resulting in a series 
of electron captures, neutron emission, and pycnonuclear reactions, through the depth of the crust, 
until the matter reaches the density of undifferentiated equilibrium nuclear matter, at which point 
it is no longer constituted of differentiable nuclei. In doing so, this reaction chain deposits 1.9 MeV 
per nucleon, distributed throughout the crust (HZ08). Although the amount of heat deposited into 
the outer crust dep ends on the composition of nuclei produced by the fusing of accreted H and He 
( Gupta et al.ll2007l ). the total amount of heat is dominated by neutron emissions and pycnonuclear 



reactions in the inner crust, giving a total heat deposition that is relatively insensitive to the nuclear 
history of the accreted matter (HZ08). 
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In a steady-state, these reactions give rise to a time-average luminosity, which is directly 
proportional to the time-average mass accretion rate (BBR98): 



L = 9xl0 32 



(M) 



Q 



10- 11 M e yr" 1 1.5MeV/amu 



erg s 



(1) 



where Q is the average heat deposited in the NS crust per accreted nucleon. Of co urse, nuclear 



burn i ng in the NS atmosphere wil l produce not pure iron, but a wide range of nuclei (jSchatz et al 



1999 



2001 



); this has been shown (jHaensel Zdunikll2003l ) to possibly change Q to range between 
1.1-1.5 MeV per accreted nucl eon. Moreover, the uncertain time-scale for the pycnonuclear burning 
of 34 Ne (lYakovlev et alJl2006l ). which arises in the Fe-chain, leaves uncertain whether that reaction 
is followed, or some other reaction takes its place on a shorter time-scale, at slightly greater density. 
However, in the final analysis, the amount of heat Q deposited in the crust is roughly equal to the 
average mass defect of nuclei at the top of the crust (near densities ~ 10 9 gcm~ 3 ) with undiffer- 
entiated equilibrium nuclear matter at the bottom of the crust, near nuclear saturation density 
2.6xl0 14 gcm~ 3 (HZ08). 

Some qLMXBs exhibit sig nificantly less photon luminosity than predicted byEq. (TT]). In severa l 



cases, such as KS 1731—260 (IWiinands et al. 



2001 



Rutledee et al.1 l2002bl : ICackett et all 12009 ) 



and 1H 1905+000 (jjonker et al.ll2007bl ) the long-term mass-transfer rate is not well constrained. 
For the accreting millisecond pulsar SAX J1808.4— 3657, however, predictions of (M) from the 
integra ted fluence agree with predict ions based on the orbital period and gravitational radiation 



losses ( Bildsten h Chakrabartv 



2001 



For this system, the low quiescent luminosity is consistent 
with strong neutrin o emission from the core, such as from nucleon or hyperon direct Urea processes 
(|Heinke et al J 120071 ). 



The energy source of the quiescent thermal luminosity aside, the therma l spectrum can be 
explained (BBR98) as due t o a realistic neutron star hydrogen atmosphere model ( Zavlin et alll996 ; 



Raiagopal &: Romani 



1996 



Heinke et al.l 120061 ) . rather than a black-body approximation imposed 



in previous works. In the case of qLMXBs, the atmosphere has been accreted from the low- 
mass post-main sequence stellar companion (B BR98) . When accre tion shuts off, metals will settle 
gravitationally on a time-scale of ~ seconds (|Bildsten et al. 1992 ). leaving a pure H-atmosphere. 
The emergent photon spectra of H atmosphere NSs is physically d ifferent - although parametrically 
similar - to that of blackbody thermal spectra (jZavlin et al.l ll996). which dramatically changes the 
derived emission area radii from the ^lkm derived fro m the black-body assu mption, to radii which 
are consistent with the entire area of the neutron star (jRutledge et al.lll999l ). 



The DCH luminosity/H atmosphere spectral interpretation has been applied in the detec- 
tion (and non-detections) of a large number of historically transient LMXBs, in cluding: Cen 
X-4 rtCanrnana et alJ l200cL iRutledee et al.1 hoOlbl): Aa\ X-l jRutledge et al.l boOlal ) 4U 1608-522 



Rutledee et al 



in't 



Zand et al 



199flh : 4U 2129+47 feutledge et akhoool): the trans ient in NGC 6440, S AX J1748-202 1 



200ll ): X1732-304 in Terzan 1 jwiinands et alU zQoj); XTE J2123- 058 fabmsick et al 



20041 ): EXO 1747-214 (jTomsick et all 120051 ): MXB 1659-29 (ICackett et al.1 12009 ): 1M 1716-315 
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(|jonker et al.l l2007ah: 2S 1803-24 5 (jCornelisse et all 12007V ): 4U 1730-22 (jTomsick et alJ 120071 ): 
1H 1905+000 (jjonker et al"1l2007bh . 



The spectral analysis of qLMXBs with the H atmosphere interpretation leads to the determi- 
nation of -Rqo, the emission radius (also called, the projected radius) of the neutron star. It differs 
from the physical radius by a factor g r , the gravitational redshift: 



-Rns 



#NS 1 



2GM NS y 1/2 
Rns c 2 J 



(2) 



The effective temperature is also affected by the gravitational redshift through the relation = 
g r T c f[. A major uncertainty in accurate measurements of Roq remains the distance D to the 
system. While high signal-toi-noise (S/N) X-ray spectra are capable of determining R^/D to 
statistical certainty, the 10-50 per cent systematic uncertainties in the distances to field sources 
directly contribute an additional 10-50 per cent uncertainty in R^. The observational motivation 
for measuring R^ - placing observational constraints on the de nse matter equation of st ate - 
requires measurements of Roo to ~ 5 per cent accuracy, or better (jLattimer &: Prakashl 12004 ). 



The observational solution is to discover new qLMXBs, the distances to which are known 
with greater certainty than the present group of field sources. Globular clusters (GCs) are the 
obvious place to search for these, as they host an over-abundance of LMXBs and their distances 
are precisely known or measurable. Historically, it was suggested that X-ray binary systems i n GCs 
are fo rmed by capture from the remnants of massive stars, i.e. neutron stars or black holes (I Clark 
19751 ). It was also noted that the ratio of X-ray sources to the surrounding stellar density is two 
orders of magnitude larger for GCs than for the Milky Way. For example, 12 out of ~ 100 known 
bright X-ray sources were lying w ithin GCs while the population of GC represents only 10~ 4 of 
the total mass of the Milky Way ([Hut et al.lll992l ). The encounter rate is related to the physical 
properties of GCs, so that the high s tellar density in their core creates a propitious environment for 
large encounter rates (|Verbuntll2002l ) . Recently, a correlation was sh own between the GC encounter 
rate and the number of hosted X-ray sources (jPooley et al.l 120031) or between the GC encounter 
rate and the number of hosted accreting NS binary systems (jHeinke et al.1 l2003d). In addition - 
cluster metallicity can play a significant role in the numb er of LMXBs in GCs (jKundu et al.ll2002l ; 



Maccarone et al. 



2003 



Ivanova 



2006 



Sivakoff et alJl2007l ). 



To date, only a handful of qLMXBs in GCs c apable of producing high S/N X-ray spectra 



are known: in Omega Ce n (jRutledge et al 



(Hein 



ce et al 



M28 ( 



2003b, 



Becker et al.l 



200S 



2002a 



2006); in NGC 6397 (Grindl av et al.1 2001bT ): in M13 dGendre et alJl2003ah : in 



Gendre et al.1 |2003bJ) ; X5 and X7 in 47 Tuc 



1); A-l in M 30 (|Lugger et alJl2007i ): in M80 (|Heinke et all l2003al ) : in M55 



and NGC 3201 (IWebb et al.1 120061 ) . Others have been reported in the literature, but with lower 
S/N, or high values of Nh, which makes them significantly less useful for the primary motivation 
for searching for qLMXBs: providing measurements of i?oo. 

In each case, the qLMXB was identified by its X-ray spectrum alone, consistent with a H- 
atmosphere neutron star with radius of ~ 10 km, without observation of an optical/IR counterpart, 
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or evidence o f an historical X-ray outburst. In two cases - X5 in 47 Tuc ([Edmonds et al.ll2002l ) and 
Omega Cen ([Haggard et al.l 12004 ) - the optical counterpart of the binary was identified following 
the spectral classification, lending support to the assertion that X-ray spectral classification can 
identify qLMXBs. 

Some qLMXBs also exhibit a power-law (PL) spectral component which dominates at photon 
energies above 2keV. Th e first detection of such a component was made on the qLMXB Cen X-4 
(jAsai et al.l Il996bl . [1998), where the power-law contributes about 40 per cent of the 0.5-10 keV 
flux. Observations of 47 T uc (qLMXBs X5 and X 7) also showed weak power-law components (with 
photon index a = 2.6—3) (jGrindlav et al.ll2001al ). However, analysis o f more sensitive observations 
concluded that PL components are not required (jHeinke et al.ll2003bl ). A more recently identified 
candidate qLMXB in 47 Tuc, has s hown a power-law component that represents 65^2 P er cen t 
of the 0.5-10 keV unabsorbed flux (jHeinke et alJbood ). NGC 6440 have possible qLMXBs with 
significant PL components, but the large absorption in the direction of this GC prevents a justified 
classi fication of the X-ray s ources due to the low photon flux where the thermal spectrum is domi- 
nant 



in't Zand et al 



2001 



). In general, qLMXBs in GCs display no or little power-law components 
while qLMXBs in the field more frequently tend to have a significant power-law component. A pro- 
posed interpretation r elates the power-law emission as the residual of a recent accretion episode 



(jGrindlav et al.ll2001al). An alternate explan ation evokes it to shock emissions via the emergence 



of a magnetic field (jCampana fc Stellal l2000i ) . However, neither of these two hypothesis are yet 
supported by observational evidence. 

A program of short-integration XMM-Newton observations, using the pn camera, have been 
undertaken to survey GCs for qLMXB candidates, to increase the number of known qLMXBs in 
GCs. This paper reports the discovery of three candidate qLMXBs in NGC 6304 from spectral 
classification, XMMU 171433-292747, XMMU 171411-293159 and the low signal-to noise source, 
XMMU 171421—292917. These were found in short XMM-Newton /pn exposures targeted specif- 
ically for the discovery of examples of this class of objects. The organisation of this paper is as 
follows. In §2.11 and §2.21 the data reduction, the observations and the astrometric analysis are 
described. In §2.31 the spectral analysis of the detected sources are presented with detailed expla- 



nations for interesting sources, including a hard source with the photon index a = —2.0 



-1.2 

-2.2- 



In© 



the observational properties of the three candidate qLMXBs are described. © covers a comparison 
with a previous observation of NGC 6304 with ROSAT ( §4.ip and a comparison with the expected 
number of qNS in GCs ( §4.2p . Finally, ©summarizes the results of this observation and concludes. 



1.1. NGC 6304 

The globular cluster NGC 6304 is located at the position R.A.=17hl4m32.1s and Dec. = 
— 29°27'44.0", (J2000) approximately 2.2 kpc from the galactic centre. Its integrated visual magni- 
tude is V = 8.22 and it has a color excess E{B — V) = 0.53. NGC 6304 is moderately compact, 
with the core radius r c = 0.21', the half mass radius thm = 1-41' and the tidal radius r% = 13.25'. 
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It is a non core-collapse GC with a concentration c = log (rt/r c ) = 1.80 (assuming a King's profile) 
and with a central luminosity density log (po) = 4.41 L@ pc~ 3 . Finally, its met allicity [m/H] is 
-0.56, implying Z = 0.00488, where Z Q = 0.0177 was used (jMontalban et al.1 12004 ). All properties 



come from the updated catalog of GCs ([Harris! Il996h or from a more recent work (jValenti et al 
20071 ). A more recent measurement of the cluster distance modulus , (m — M) = 13.88 ± 0.03, 
provides an heliocentric distance of 5.97±0.08kpc (jValenti et al.l l2005). 



2. Data Reduction and Analysis 
2.1. Observation 

The target was observed on 4 Sept 2006 using XMM-Newton. The three X-ray CCD cameras 
- MOS1, MOS2, and pn - were operating in imaging mode, with exposure times of 12618 sec, 
12626 sec and 11037 sec, respectively. All exposures were m ade with the medium filter. The present 
analysis is focussed on the pn camera (jStriider et al.ll200ll ). due to its greater sensitivity in the low 
photon energy range (0.5-2.0 keV) where qLMXBs are brightest. However, a drawback of the pn 
instrument is the multiple bad columns, which are accounted for during the analysis. 

Data reduction was performed using XMM-Newton Science Analysis System v7.0.0, using 
standard analysis procedures. The command epchain performs all required steps for preliminary 
data reduction. The data were checked for proton flares, but the standard analysis found none. 
However, visual inspection of the exposure reveals a statistically significant background variability 
during the observation. Despite this, the full observed time interval has been used for this analysis. 
The spectral quality of the data is quantified using the task epatplot, as recommended before 
choosing the energy band0. Over the whole field, epatplot shows a significant discrepancy below 
0.45 keV between the distribution of the various pattern events (single and double) as a function 
of energy and the expected curves, therefore excluding the use of this low energy band. The 0.5- 
lOkeV spectral energy range was therefore chosen during the data reduction, the detection and 
the subsequent steps of the analysis. Only the event patterns recommended by the XMM-Newton 
Science team were used, that is, single and double patterns for the pn camera, and single, double 
and quadruple patterns for the MOS cameras. An exposure ma p was then created fo r source 
detection, which was performed with wavdetect from CIAO v3.4 (jFruscione et al.ll2006l ). Using 
the pn data, a total of 11 sources were detected above a minimum wavdetect relative exposure 
of 0.2 (see Table [Hand Fig. Q]). Other wavdetect parameters used are the wavelet scales (1.0 2.0 
4.0 8.0) and the significance threshold (3xl0 -6 ). Using the same detection parameters, 6 other 
sources are detected on a MOS1+MOS2 combined image, most of them being coincident with bad 
columns on the pn camera. The source ID numbers in this paper (from to 10) are given by the 
default numbering of wavdetect, attributing the ID numbers by decreasing right ascension. The 



^rom XMM-Newton Science Operations Center XMM-SQC-CAL-TN-0018 lGuainizzil (|2008h 
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additional MOS sources are labelled Ml to M6. The pn data was used to address the primary 
purpose of this observation. The source detected on the pn camera were spectrally fitted in order 
to identify candidate qLMXBs. For completeness and as a secondary objective of creating a catalog 
of X-ray sources in the field of NGC 6304, the sources detected on the MOS1+MOS2 image are 
reported in Tabled 



2.2. Astrometry 

The astrometric positions given in Table [T] were derived in the following manner. Astrometry 
was first obtained using wavdetect with the pn data, the MOS1 data, and the MOS2 data, all 
separately. Comparing the relative astrometry for sources detected by each of these cameras, 
discrepancies are noticed in the relative astrometry for some sources. Sources #5, #6 and #9 
relative distances from a reference source (see below) in the pn camera, compared with that from 
the MOS cameras, were not consistent with the statistical positional uncertainties. The observed 
inconsistencies seem to be due to bad columns on the pn camera, which are not accounted for in 
the wavdetect algorithm; thus wavdetect misestimates their positions in the direction opposite to 
the bad column (which is the direction of the discrepancy in all three cases). 

Therefore, combining astrometric data from the pn and MOS1 cameras produced the X-ray 
astrometric positions of the detected sources. For most of the sources, the MOS1 data were used to 
deduce the positions of the sources. Some sources detected on pn, however, were not detected on 
the MOS1 camera (#2, #3, and #8). The astrometric positions of these sources are the pn position 
with a corrective offset calculated from the difference in the relative astrometry of the concerned 
source with source #1 on the two cameras. For example, RA^ osl = RA^ n + ARA^ ol is used to 
calculate the right ascension of source #3 as it would be on MOS1. The reference X-ray source is 
required to be detected on both cameras and not adjacent to a bad column. Only source #4 and 
#4 fulfil these conditions in the pn camera data, but source #4 was discarded since it is possible 
that it corresponds to unresolved multiple sources due to its presence in the core of NGC 6304. 

Finally, an absolute astrometric correction - into the frame of the Two Micron All Sky Survey 
(2MASS) - is applied through the association of the X-ray source XMMU 171411-293159 with its 
counterpart 2MASS 17141152-2931594 (99.988 per cent confidence, see ffT2jh shifting the X-ray 
positions by a change of -1.66 arcsec in R.A. and -1.02 arcsec in Dec. Previous to the frames 
alignment, the separation between the pn and the 2MASS sources, 1.77 arcsec is consistent with 
the absolute astrometric uncertainty of XMM-Newton (2 arcsec within loll). The sy stematic un- 



certainty on the position, after the correction is 1.2 arcsec (la) (j Jeffries et al,ll2006l ). Also, the 
statistical uncertainty from the source detection (up to ~ 1 arcsec, depending on the source; see 
Table QjJ should be added in quadrature to the systematic error to obtain the total uncertainty in 



2 from XMM-Newton Science Operations Centre XMM-SOC-CAL-TN-0018, 



Guainizzil l|200sh 
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position. 



2.3. Spectral Analysis 



The X-ray counts for spectra are extracted using a radius of 25 arcsec around each source, 
accounting for 81 per cent of the total energy from the source at 1.5 keV. The background counts 
are extracted using an annulus centred at the position of the peak of the exposure map (that is, 
the image centre), and with an average radius equal to the distance between the source being 
analysed and the image centre. An annulus width of 40 arcsec is used for all sources, except 
those sources close to the centre (sources #4 and #5) for which a larger background area (100 
arcsec circles) was extracted. Around each detected source, a circular area of 35 arcsec radius are 
excluded from the background annulli. A circular area of 35 arcsec excludes 87 per cent of the 
overlapping source counts (at 1.5 keV). In all cases, the sources contributes ^ 30 per cent of the 
Poisson uncertainty in the number of background counts - to the background regions. This choice 
of background area addresses mirror vignetting, which is the main cause of background variation 
across the field. At the source count rates encountered here, the effect of pile-up can be neglected. 
After the counts extractions, response matrices were created using rmf gen and arf gen, rather than 
using the predefined ones, in order to take into account the enclosed energies and the extraction 
radii. 

A threshold of 70 counts above the background was imposed before considering the X-ray 
source for a full spectral analysis. Out of 11 detected sources, six satisfy this cut (ID #1, 4, 5, 6, 
9 and 10) while the remaining five are too faint for a full spectral analysis (ID #0, 2, 3, 7 and 8). 
The counts were binned into energy bins between 0.5 and lOkeV with 25 counts per bins (or 20 
counts per bins for sources with less than 100 counts in excess of the background). 



Using XSPEC vl2.3 (|Arnaudll996T). the extract ed spectra are fitted with a tabulated neutron 
star hydrogen atmosphere model (jZavlin et al.lll996l ). For this model, the radius is measured from 



the normalisation parameter (Roo/ (<i/10kpc)) 2 , assuming a distance d = 5.97 kpc. If this model 
provides a statistically acceptable fit (Null hypothesis probability >10 -2 ) and values of and kT c g 
in the range of previously observed qLMXBs (see Table H]), 5 — 20 km and ~ 50 — 150 eV respectively, 
then, the source is classified as a candidate qLMXB. The spectral parameters are then confirmed 
with a singl e- comp onent mode l of neutron star hydrogen atmosphere included in XSPEC (NSA 



model: nsa, Izavlin et all (Il996h ). Alt ernatively, two other NS atmosphere models are available i n 



XSPEC: nsagrav (jZavlin et al.l 11999 ) and nsatmos (IMcClintock et al.1 12004 : Irleinke et al.1 12009 ) . 
When using any of those models, the mass is kept constant at 1.4M and the model normalisation 
parameter (l/d 2 ) is set assuming d = 5.97 kpc. The nsa model contains a magnetic field parameter 
that is set to zero here while nsagrav and nsatmos assume non-magnetic atmospheres. In all cases, 
the models nsa, nsagrav, and nsatmos give consistent parameters. Regarding the surface gravity 
g, nsa has been implemented with a fixed value g = 2.43x10 cms" 2 , nsagrav has been computed 
for a range g = (0.1 — 10) x 10 14 cms -2 and nsatmos in the range g = (0.63 — 6.3) xlO 14 cms -2 
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In the case where the NSA tabulated model provides a non-acceptable fit, a visual inspection 
and a F-test indicate whether a photon power-law should be added to the spectral fit. If a significant 
portion of the high energy tail is in excess of the NS atmosphere fit, a simple power-law is added. 
This is justified by the fact that qLMXBs sometimes display a hard power-law component which 
dominates the spectrum above 2keV. If the fit is still not acceptable, the spectrum is fitted 
with a power-law alone and the flux is obtained. For low signal-to-noise ratio (S/N) source, the 
slope of the power-law is fixed at a = 1 to estimate the flux. For all model fits, the photo- 
electric galactic absorption (the multiplicative model wabs in XSPEC) with a fixed value of Nh_ = 
0.266 xlO 22 atoms cm (written Nh22 = 0.266 afterwards) at the position of the NGC 6304_l is 
included in the spectral fit. 

A major technical problem arises during the spectral fitting with the nsa model, as well as with 
the nsagrav and nsatmos models. All three of them were implemented in XSPEC with restrictive 
parameter spaces. In particular, the radius parameter space is defined only between 5 km and 
20 km for nsa, between 6km and 45 km for nsatmos and between 6 km and 30 km for nsatmos. 
In this case, while the best-fit may be statistically acceptable, and the parameter values obtained 
are consistent with those expected for a qLMXB, the parameter uncertainty range is greater than 
that tabulated in the XSPEC models and the uncertainty calculation fails due to zero diagonal 
components in the parameters matrix. The consequence of the imposed cut off in the parameter 
space would be a misrepresentation of the parameter uncertainties, pr eventing any concl usions to 



be drawn from the spectral fit. In those cases, the tabulated model (jZavlin et al.lll996l ) and the 
parameters resulting from the fit are used, since they do not suffer from this constraint. In all cases 
of identified qLMXBs, the MOS1 and MOS2 spectra are extracted in a similar manner as described 
above, and used those to perform a simultaneous fit. This procedure reduced the uncertainties on 
the temperature and radii of the candidate qLMXBs, by up to fifty per cent. 

The results of this spectral analysis are given in Table [3j The following subsections present 
two sources that are not candidates qLMXBs but which were above the signal-to-noise cut-off for 
spectral analysis and the results of the spectral fitting of source #8, an X-ray source with an 
extremely hard photon power-law. Finally, the spectral results for the three qLMXB candidates 
are described in 331 



2.3.1. XMMU 171516-292224 - Source #0 

This source, located far from the optical centre of the GC (~ 52r c ), is not spectrally consistent 
with a qLMXB at the distance of NGC 6304. Using the tabulated H-atmosphere model with a fixed 
Nh(Nh,22 = 0.266), gives a non acceptable fit (x^/dof (prob.) = 3.0/12 (3xl0 -4 )) to the spectrum 
of the source. Allowing Nh to vary provides an acceptable fit x^/dof (prob.) = 1.2/11 (0.28) with 



3 From 



http://cxc.harvard.edu/toolkit/colden.jsp 



using the NRAO data 



Dickev fc Lockmanl [|l990h 
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Nh, 22 = but the resulting projected radius Roc = 1.3+Q22 ^ m excludes a possible quiescent 

neutron star interpretation since the value is significantly smaller than previously measured for such 
objects. 

A power law spectral model is therefore chosen to characterize the X-ray emission and it is 
found that with Nh, 22 = 0.266, the photon index is a = — 0.24^gg 2 for which the statistic is %^/dof 
(prob.) = 1.0/12 (0.45). If the absorption is chosen to be a free parameter, the fit is still statistically 
acceptable (xl/doi (prob.) = 0.84/11 (0.60)) with the parameters N H ,22 = and a = 0.3^- 

Finally, a search for a possible 2MASS counterpart found no source within 3 sigma of the position 
of source #0. 



2.3.2. XMMU 171453-292436 - Source #1 

Using the H-atmosphere tabulated model, a non acceptable fit with x^/dof (prob.) = 2.47/15 
(0.001) is obtained when Nh is fixed at Nh,22 = 0.266. On the other hand, as Nh is allowed to 
vary, the best fit (x^/dof (prob.) = 0.98/14 (0.47)) projected radius = 0.281^8 km excludes 
the qLMXB classification for this source. In addition, the best-fit (with fixed Nh) exhibits an 
excess high energy tail, suggesting an additional power-law component at high energies. Adding 
a power-law component to the assumed spectrum results in a dominating power-law (~ 100 per 
cent of the source 0.5-10 keV flux) of slope a ~ 0.5. The NS atmosphere component is therefore 
abandoned and the spectral fit is restricted to a single power-law for which the best-fit photon 
index is a = 0.52^Qg 5 with a value x^/dof (prob.) = 0.96/15 (0.5) in the case where Nh,22 = 0.266. 
However, as Nh,22 = 2.3^}'g, the photon index is a = l-7lg 9 f° r w hich x^/dof (prob.) = 0.65/14 
(0.82). 



2.3.3. XMMU 111411-293222 - Source #6 

The NSA tabulated model produces an acceptable fit (x^/dof (prob.) = 1.3/12 (0.21) for 
N H ,22 = 0.266 and x^/dof (prob.) = 1.3/11 (0.21) for N H ,22 = 20.8l^ 8 8 ) but in both cases the 
obtained parameters (Roo < 1km and kT e Q ~ 0.5 keV) are not in accordance with expected values 
for NSs in quiescence. Also, visual inspection of the fitted spectrum indicates a significant high 
energy tail. The fit of the model (tabulated model+PL) is also statistically acceptable (x^/dof 
(prob.) = 0.96/10 (0.48)) but, similarly to source #1, ~ 100 per cent of the unabsorbed flux comes 
from the power-law component. The single power-law model (a = 0.64to'go) provide a statistically 
acceptable fit (x^/dof (prob.) = 0.75/12 (0.70)). The spectrum fit with a fixed slope (a = 1) 
results in an almost unchanged Xu statistic (x^/dof (prob.) = 0.75/13 (0.71)). Finally, as both the 
absorption and the photon index are allowed to vary, the best fit parameters are a = 0.45l£;jj and 
N H ,22 < 0.9 for x^/dof (prob.) = 0.78/11 (0.66). 
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2.3.4- XMMU 171413-293415 - Source #8 

Even though it is one of the faint sources, XMMU 171413—293415 is interesting since a single 
power-law spectral model suggests a spectrally hard source. As Njj is kept fixed at Nh,22 = 
0.266, the best fit photon index a = -2.0+2:2 with X^/dof (prob.) = 1.726/7 (0.10) implying 
an unabsorbed flux Fx = 2.8xl0~ 13 erg cm -2 s _1 . When the absorption is a free parameter, the 
photon index, a = —2.1, is similar to the case where Njj is fixed , but the large xt value, x^/dof 
(prob.) = 2.02/6 (0.06), greater than 2 does not allow for error estimation. 

Most of the counts are in the range 5-10 keV (Fig. [2]). The 5-10 keV light curve (100 sec bins) 
exhibits no variability (Fig. [3]). The average number of counts per bin is 1.12, and 0.78 counts are 
expected on average due to the background. So, the average number of source counts is 0.43. A 
bin with 5 counts is found and the Poisson probability of finding a bin with such a peak in the 
lightcurve with a mean of 1.21 is 2.2 per cent. Also, the peak does not represent a variation by more 
that a factor of 10 of the mean number of source counts. The spectrally hard source is consistent 
with constant emission throughout the duration of the observation, to within a factor of < 10 in 
any 100 second bin. This excludes a source type such as an X-ray flash or soft gamma-ray repeater 
(SGR), which tend to be spectrally hard in the 1-10 keV energy band as exhibited by this object, 
but exhibit transient emission of duration ~few msec, to ~100 sec. No detailed variability analysis 
can be performed due to the highly variable background. In conclusion, no variability is detected 
on a < 100 sec time-scale. 

To confirm the existence this source, the MOS2 camera data were used (the CCD 6 on MOS1 
containing source #8 was inoperative during the observation). The same data reduction as de- 
scribed for the pn camera and the same detection procedure is performed. As explained before, 
the data quality of the MOS camera is much better than the pn but its sensitivity is not as good 
as the pn camera. On MOS2, a total of 73 counts are detected within 25 arcsec of the expected 
position of source #8, with 2254 ± 47.5 counts in the background annulus. From the scaling 
factor between the source and the background areas, 50.3 ± 1.1 background counts are expected 
within the source extraction region. Therefore, there is an excess of 22.7 ± 8.6 counts, corre- 
sponding to a low-significance detection (2.6<r); the corresponding flux, assuming a power-law of 
slope a = —2.0, is 2.6xl0 _13 erg cm -2 s _1 using the tool webPIMMSn, which is consistent with 
the 2.8xl0 _13 erg cm" 2 s _1 value measured from the pn spectrum. Cross identification with the 
2MASS and Digital Sky Survey (DSS) catalogues did not exhibit any counterparts in the J and B 
bands within 3<r of the X-ray position, respectively. 

While the low signal-to-noise ratio of this detection (S/N = 4.0) does not permit any detailed 
analysis of this source, it is unusually spectrally hard for an object detected in the 0.5-10 keV band. 
Examination of the unfiltered event list (i.e. including all event patterns) finds no unusual number 
of rejected events, which does not support an interpretation of the source detection as due to a 



4 available at |http://heasarc. nasa.gov/T ools/w3pimms. html 
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particle impact on the detector or an instrumental effect such as a transient thermal excess on the 
detector. Higher sensitivity observations are required to better characterize this source. 

3. Candidate qLMXBs 

This section presents the spectral analysis of sources which, on the basis of this spectral 
analysis, are candidate qLMXBs. 

3.1. XMMU 171433-292747 - Source #4 

This X-ray source, located near the centre of the cluster, has the highest count rate among the 
sources in the field (43.8 ± 2.8ctsks~ 1 ). 

The fit with the tabulated NS atmosphere model is statistically acceptable (x^/dof (prob.) 
= 1.24/21 (0.20)) and is producing values for i?^ and kT e s in the range expected for a NS in 
quiescence. However, there exists a systematic excess of counts above the best-fit model at photon 
energies larger than 3keV suggesting an additional high-energy spectral component, which are 
accounted for by the addition of a power-law. The low probability value of the F-test (prob=0.014) 
suggests that adding a power-law component better describes the data than the absorbed H- 
atmosphere alone. The best-fit (tabulated) NSA model with the power-law component has a fit 
statistic x^/dof (prob.) = 0.69/19 (0.83) (see Fig.©. Also, = 8.1±|J km and kT eS = 127±|g eV 
with Nh ,22 = 0.266. These values are in agreement with other known qLMXBs (Table UJ). However, 
when leaving the absorption free Nh,22 = 0.51^q3q for this model, the fit is acceptable (x^/dof 
(prob.) = 0.63/18 (0.88)) but the radius is larger than expected for a neutron star, = 47^3,1° km. 

The fit with fixed Nh is confirmed with a nsa+powerlaw model, for which the Xu statistic 
is x^/dof (prob.) = 0.83/19 (0.67). The radius (i?NS = 8.3^3 q km) and temperature (kT e g = 
117^4 eV) suggests a candidate qLMXB since the X-ray spectrum is consistent with the thermal 
spectrum of observed qLMXBs at the distance of NGC 6304. The slope of the power-law is a = 
1-5^08 an d this component of the model contributes to 49 per cent of the total unabsorbed 0.5- 
lOkeV flux of the source. A simultaneous fit using the MOS1, MOS2 and pn spectra improved the 
statistics (x^/dof (prob.) = 0.85/42 (0.75)) and the uncertainties on the parameters. The radius 
and temperature of this candidate qLMXB are now 

Rns = 8.1 km and kT eS = 122±jg eV, and 
the photon index is a = 1.2+q' 8 . The value i?NS is converted to Roo assuming a 1.4 M neutron 
star using Eq. [2j R^ = 11.6±£jj km. 

The absence of significant intensity variability over the time scale of the observation further 
supports the classification of this source as a qLMXB. Moreover, the luminosity is consistent with 
being the same as that observed during the ROSAT observation 14 yr earlier (see Table [5] and 
section |4~T1) . This is consistent with the expected stable thermal luminosity from a qLMXB on this 
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time-scale fBBR98. lUshomirskv k RutlededfeOOlh . 

The possibility of two unresolved sources - one spectrally hard and one spectrally soft - 
is investigated by separating the data into soft (< 1.5 keV) and hard (> 1.5 keV) images and 
performing a source detection as described in Sec. 12.11 Source #4 is detected in both bands, 24.6a 
and 4.5<7 significances in the soft and hard bands respectively (the 1.5 keV division was selected 
as the lowest possible photon energy which provides a significant source detection - a cutoff of 
2.0 keV finds no source with > 3a significance above 2.0 keV). A positional offset of 3 ±0.9 arcsec is 
measured between the positions of the source in the two images, which is marginally consistent with 
a single source. It is possible that an unresolved hard X-ray source lie close to the candidate qLMXB 
and was interpreted as a high energy tail in the spectrum of source #4. In other words, it is not 
possible to exclude that higher spatial resolution observations may resolve XMMU 171433—292747 
into multiple X-ray sources. 



3.2. XMMU 171411-293159 - Source #9 

3.2.1. Spectral Characterisation of the X-ray source 

The fit with the tabulated NS atmosphere model with fixed Njjis statistically acceptable and 
does not require an additional power-law (Fig. EJ). The xl statistic is x^/dof (prob.) = 1.29/16 
(0.20), the projected radius is Rqo = 15.3ti; 5 2 5 km and the effective temperature is kT e ^ = lOOl^g eV. 
Letting the absorption vary also gives a statistically acceptable fit (x^/dof (prob.) = 1.38/13 
(0.16)) with the following parameters: N Hj2 2 = 0.20l£f£, = 9.8±$g km and kT cS = 115±§J eV, 
all statistically consistent with the fit with a fixed absorption. 

The confirmation with any of the models nsa , nsagrav and nsatmos fails because the fit, while 
statistically acceptable, produces error regions larger than the allowed parameter space in XSPEC. It 
is therefore not possible to derive the uncertainty region for the spectrum of XMMU 171411—293159 
using these models. In consequence, the results of the tabulated model show that the model fit is 
consistent with the thermal spectrum of a qLMXB at the distance of the NGC 6304. Moreover, 
no significant variability over the time-scale of the observation is observed, with weak constraints. 
To better improve the statistics and uncertainties, the MOS1, MOS2 and pn spectra are fitted 
simultaneously. The resulting fit is statistically acceptable (%^/dof (prob.) = 1.20/31 (0.21)) 
and the obtained parameters are: kT e g = 115^^ eV and Roo = 10.7i 3 '?km. Those values for 
XMMU 171411-293159 are reported as a candidate qLMXB in Tables [3] and H 

For completeness, additional spectral fits are provided. First, A single power law with fixed 
absorption {Nh,22= 0.266) gives an acceptable fit (x^/dof (prob.) = 1.37/14 (0.16)) with a soft 
photon index, a = 3.511"q' 3 °. Leaving the galactic absorption Nh free Nh,22 = 0.57^q'^ 7 , a steeper 
power-law is found a = 5.0j^ 5 with a x^-statistic Xv/dof (prob.) = 1.28/13 (0.21). A single 
temperature Raymond-smith plasma fit with variable absorption and solar metallicity also produces 



-14- 



an acceptable fit (x^/dof (prob.) = 1.03/13 (0.42)) with the following best-fit parameters: Nh,22 = 
0.19±£^ and kT cS = 0-76±g;f| keV. Thus, the spectrum of this (highest S/N) X-ray source is 
consistent with other spectral interpretations; however, such a steep power-law (either a = 3.5, or 
a = 5.0) is usually interpreted as indicating a thermal spectrum. 

As mentioned above, the pn camera suffers bad pixel quality. Two of the detected sources (#6 
and #9) are overlapping with bad columns. So for XMMU 171411—293159 , an additional analysis 
using the MOS 1 and 2 cameras is performed. Source #9 is detected with a signal-to-noise ratio 
of 17.7 and 17.6 for MOS1 and MOS2 respectively. Using the tabulated model, the spectral fits of 
the MOS1 and MOS2 data and the simultaneous spectral fit of the MOS1, MOS2 and pn data are 
all statistically acceptable and give results similar to the ones obtained from the pn camera alone. 



3.2.2. Search for an IR counterpart 



Using the online 2MASS Point Source Catalogue (2MASS-PSC), a probable IR counterpart 
(2MASS 17141152-2931594; mj = 8.796(22), m H = 8.361(18), m K = 8.213(20), where the num- 
bers in parenthesis are the la uncertainties in the preceding digits) is found located at a distance 
of 1.77 arcsec (before astrometric correction to the 2MASS frame). For comparison, there are 6 
stars of this magnitude or brighter in an annulus (r; n = 4.5' and r out = 8') about the centre of 
the GC. The probability that another source as bright or brighter lies as close or closer to our 
X-ray source is 0.012 per cent. Therefore, the association between XMMU 171411—293159 and 
2MASS 17141152-2931594 is identified on t he basis of spatial proximity, with 99.988 per cent 
confidence. Also, the USNO-B1.0 catalogue ( Monet et alJl2003l ) lists an object at this location 
(separated by <0.5 arcsec) with B = 11.66 and B = 11.28. The two different magnitudes may 
indicate time variability of the counterpart. 

A photometric study of NGC 6304 finds an observed V = 10.65 for 2MASS 17141152-2931594 
(S. Ortolani, private c ommunication, 2008) ; this is about 3 magnitudes above the observed tip of 
the red giant branch (jOrtolani et al.l |2000| ) (see Fig[7|), however, due to saturation of the CCD, 
the photometry of very bright stars, like this counterpart, was not performed - thus, objects as 
bright as this were not included in the published dat4j- Inspecting by eye a U-band exposure of 
NGC 6304 (S. Ortolani, private communication, 2008) shows that ~ 15-20 such saturated stars were 
not included in the BV photometric study of this GC. 

The color-magnitude diagrams (Figs. [6] & [7]) contai ning theoretical isq chrones for ages in the 
range t = 10 8 - 55 - 10 iai5 yr and metallicity Z=0.00488 (|Marigo et alJbood ) show that the optical 
counterpart is consistent with a post-asymptotic giant branch (post-AGB) star, both in V and J 
bands. Even if the time scale of this stage is short, finding such an object in a GC is not completely 



5 The photometry for 2MASS 17141152—2931594 was performed by S. Ortolani from a separate shorter (10 sec) 
V-band exposure 
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unexpected. Based on pre vious work, 16 plan etary nebulae (PN) are expected to be found in the 
GCs system of the galaxy (jJacobv et al.lll997l ). This corresponds to 6.7x 10~ 7 PN L _1 , using the 
total lumino sity of the GCs sy stem, ~ 2.4x10 Lq (SeckerJ [l992j). Using the estimated mass of 
NGC 6304 (|Gnedin et alJbood ) @ and the average mass-to-light ratio of GCs, M/L = 1.7 dCaputo 
1985), 0.25 PNs are expected in this cluster. Since the lifetime of a PN, ~ 10 yr ( Jacobv et al 



1997 ). is comparable to that of a post-AGB object, ~ 10 3 — 10 4 yr ( Siodmiak et al. 20081 : 



Blocckcr 



19951 ). finding one post-AGB star in NGC 6304 is not completely unexpected, and is consistent 
with the object classification. 

To obtain the bolometric luminosity of 2MASS 17141152-2931594 at the dist ance of NGC 6304 
the r eddening E(K - V) is first calculated using E(K - V) /E(B - V) = -2.744 rtllieke & Lebofskv 
198fil ). where E(B-V) = 0.53 for NGC 6304. The intrinsic color is (V-K) = ( V-K)+E(K-V) = 
2.38 — 1.45. Using bolometric corrections obtained from 2MASS photometry (IMasana et al.1 12 006) 
as a f unction of (V — K)q and of the average metallicity of the cluster [m/H] = —0.56 (jValenti et al. 
20071 ) . the i^-band b olometric correcti on is BCk = 0.892. Using a i'T-band extinction correction, 
A K = 0.36 E(B-V) (|Fitzpatricklll999l ). m hoi = m K -A K + BC K = 8.914 is finally calculated to be 



able to work out the bolometric luminosity. Taking a z ero absolute b olometric magnitude Mbd = 
star to correspond to a luminosity of 2.97x 10 35 erg s _1 ( Harwitll2006l ) and the distance modulus, the 
absolute bolometric magnitude is M^ \ = —4.96, for a bolometric luminosity of 2.9xl0 37 erg s _1 . 

As a check on the value calculated for the bolometric luminosity, a less precise method of 
estimation is used (since it does not take into account effects of realistic atmospheres, such as limb 
darkening), by fitting a blackbody curve to the unabsorbed photometry in BVJHK bands. For a 
blackbody source at the distance of NGC 6304, this resulted in a best fit (x 2 /dof = 5.8, for 3 dof) 
of R = 46±l.Oi? , T cS = 6620±140K, for a total luminosity of L = A-kR 2 oT^ = 1.4xl0 37 erg s^ 1 ; 
this is a factor of x2.1 below that derived through the more precise method, which may be indicative 
of the systematic uncertainties in approximating the total luminosity using a blackbody fit to the 
BVHJK data, but which is regarded as consistent with the more precise BCk calculation method. 



3.2.3. Examining the Hypothesis that XMMU 171411-293159 is a Foreground Star 

The hypothesis that 2MASS 17141152-2931594 is a field dwarf star is investigated. First, this 
star (TYC 6824-713) is listed in the Tycho-2 catalog (issued from a re-analysis of the Hipparcos 
data) providing proper motion, but no parallax information (?). The proper motion observed for 
this star is -6.6±2.8 milli-arcseconds per year, marginally consistent (at 2.4<r) with a zero proper 
motion, as would be expected for a cluster member. This does not resolve the possible cluster 
membership versus the foreground star hypothesis. 

If reddening dust is distributed uniformly between the observer and NGC 6304, the source, 



Online data at |http: / /www.astro.lsa.umich.edu7^ognedin/gc/vesc.dat 
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assuming it is a main-sequence star, is most consistent with an My = 4.4 object at a distance 
170 pc, with E(B — V) = 0.03, based on t he absolute magnitude My and color B — V values of 



stars in the Hyades (jde Bruijne et all 12001 



For main-sequence stars, the typical X-ray to bolometric flux ratio seems to reach an upper limit 
at 10" 3 jvilhu fc Walteilll987l ). For the association XMMU 171411-293159/2MASS 17141152-2931594, 
the X-ray to bolometric flux ratio is of the order of ~ 10~ 5 , therefore consistent with a main- 
sequence foreground star. 

In a B - V vs. J - K color-color diagram (Fig. EJ), XMMU 171411-293159 lies marginally 
off the main-sequence (by between 2 and 3 a) for a wide range of metallicities (Z = 0.0001 — 
0.03), lying closest to stars in the mass r ange 0.7-0.9 M@, depending on metallicity, according to 
theoretical isochrones (IMarigo et al.ll2008l ). This marginal result from the color-color diagram does 
not definitively exclude the possibility that the star is in the foreground. 



The X-ray emission from dMe stars is spectrally comparable to those of RS CVn ( Si ngh et al 



1996). Typically, RS CVn exhibit a two-temperature Raymond-Smith plasma model (jDempsey et al 
1993). Since the number of counts are not sufficient for such a fit with all parameters left free, 
the spectrum of XMMU 171411—293159 is fitted twice with a two-temperature Raymond-Smith 
plasma model keeping the high temperature component fixe d at the lowest (kT? = 0.93 keV) and 
highest value {kT 2 = 3.45 keV) from the sample of RS CVn (jDempsey et al.lll993l ). The metallic- 
ities are also held fixed at the solar value, Z = 0.0177. The normalization of the Raymond-Smith 
model in XSPEC model provides the emission measure (EM = J n e njidV), a distance dependent 
value. The best fit parameters for those two statistically acceptable fits are Njj22 
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sponding to kT2 = 0.93 keV and Nh,22 
and EM 2 /EMi < 0.55 for kT 2 = 3.45 keV. 



cm" 3 and EM 2 /EMi = ll^xlO" 5 for the first fit corre- 



0.14i°;^, kT x = 0.84l^| keV, EMi 



-0.25 ■ 



1.3±g;|xl0 56 cm" 



The value of EMi lies between 3 and 7 orders of magnitude above the typical values (EMi = 
0.1 — 3xl0 53 cm~ 3 ), for the fits performed with the maximum and minimum kT 2 respectively, 
assuming 2MASS 17141152-2931594 lies in NGC 6304. Thus, an object with a typical EMi would 
have to be at a distance less than ~ 330 pc. 

The ratio of the two emission measures EM ? /EMi , however, is a distance independent value, 
that is in all cases > 1 for coronally active stars (jDempsey et al.lll993l ). In both cases, the derived 
90 per cent confidence upper limits given before lie below values observed from coronally active 
systems. Therefore, the observed X-ray spectrum is inconsistent with the typical two-temperature 
plasma of a coronally active star, independently of the distance of the IR counterpart. 

In conclusion, the optical/IR colors and proper motion of the counterpart are in marginal 
agreement with those of a main sequence object located in the foreground of NGC 6304. However, 
the X-ray spectrum of XMMU 171411—293159 is inconsistent with those of typical coronally ac- 
tive stars. This was shown using a distance independent value. In addition, the possible optical 
variability of the star (in the 5-band, see Sec. I3.2.2|) presents further support that the star is not 



-17- 



merely a typical dwarf main sequence star. At the same time, the X-ray spectrum is consistent 
with a H atmosphere qLMXB at the distance of NGC 6304; and the V vs. B-V and J vs. J-K 
CMDs support interpretation of the 2MASS counterpart as a post-AGB star in NGC 6304, which 
is evolutionarily consistent with the X-ray spectrum of a qLMXB. Without excluding classification 
as a foreground star, we identify XMMU 171411-293159 as a candidate qLMXB in NGC 6304; 
however, an optical/IR spectrum of the 2MASS counterpart should be obtained to confirm that 
the IR source is an actual post-AGB star in NGC 6304, and not a marginally spectrally unusual, 
optically variable main-sequence star (0.7-0.9 Mq) at a distance of ~few 100 pc. 



3.3. XMMU 171421-292917 - Source #5 

The NS H-atmosphere tabulated model applies for this low signal-to-noise source. The fit is 
statistically acceptable (x^/dof (prob.) = 1.09/16 (0.36)) and the parameters are consistent with 
expected values (kT eS = 70±H eV and = 23±?g km with N H ,22 = 0.266; See FigEJ. When 
leaving N H free, the fit results become x£/dof (prob.) = 1.12/15 (0.33) with N H<22 = 0.46±°;^, 
kT eS = 49l^eV and > 8.9 km (the upper bound error estimate failed to converge), i.e. 
consistent with the previous fit when Nh was fixed. 

The nsagrav model is then used to confirm the first fit with the NS H-atmosphere model. Even 
though acceptable (%^/dof (prob.) = 1.23/16 (0.23)), the fit of this low signal-to-noise spectrum 
does not allow for representative error estimates, as explained in H2.31 The other two XSPEC neutron 
star atmosphere models also have the same issue. Also, the simultaneous fitting of the pn, MOS1 
and MOS2 spectra did not provide an improvement of the uncertainties. 

Using the online 2MASS-PSC catalogue, a possible counterpart, 2MASS 17142095-2929163 
(raj = 13.063(29), m# = 12.645(33) and mx = 12.510(35)), is identified at an angular distance of 
0.45 arcsec. Using the same calculations as for source #9, the probability of finding a counterpart as 
bright or brighter, as close or closer is 0.084 per cent when 179 stars with magnitude mj < 13.063 
are found in an annulus of inner radius 2 arcmin and outer radius 4 arcmin around the optical 
centre of NGC 63049, implying an association with 99.916 per cent confidence. 

Using the procedure described in § 13.2.21 the bolometric magnitude and luminosity of the 
IR counterpart are calculated from the V-band magnitude V = 14.85(5), the B-band magnitude 
B = 14.75(9) and the color B — V = 0.90(7) (S. Ortolani, private communication, 2008). The values 
found are mboi = f^K — A-K + BCk = 12.82, implying an absolute bolometric magnitude M^ Q \ = 
— 1.07. Therefore, the bolometric luminosity is L^ i = 7.96x 10 35 erg s _1 , when the bolometric 
luminosity for the zero absolute bolometric magnitude is 2.97x 10 35 erg s _1 . As a check on the value 



7 Prior to the astrometric correction using XMMU 171411-293159 and 2MASS 17141152-2931594, the angular 
distance between the X-ray and 2MASS sources was 1.59 arcsec. The corresponding association confidence would 
have been, in that case, 98.95 per cent. 
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calculated for the bolometric luminosity, a blackbody curve is fitted to the unabsorbed photometry 
in BVJHK bands, assuming the source at the distance of NGC 6304. This resulted in a best 



fit (x 2 /dof 
L 

calculation above. 



for 3 dof) of R = 
3.15xl0 35 erg a -1 ; 



5.87±O.14 J R , T eS = 7180±150K, for a total luminosity of 
this is a factor x2.5 below that calculated using the BCk 



If the X-ray emission is coming from the star 2MASS 17142095—2929163, the X-ray to optical 
flux ratio for this giant star is therefore Fx/F^ i = 10~ 05 . From t he ROSAT All Sky Survey cat- 
alogue of bright late- type giant and supergiant (jHunsch et al.lll998l ). the mean calculated X-ray to 
optical flux ratio for all 450 stars observed in that catalogue is Fx/ F^ = io _5 - 6±0 - 6 (la). The 0.1- 
2.4 keV range was used to estimate the X-ray flux. The 2MASS counterpart and the X-ray source 
#5 are consistent with being a giant star (2.6a). However, if X-rays from XMMU 171421—292917 
were due to a typical giant star, then similar X-ray flux from other giant stars in this GC would 
be expected; yet no other giant star on the outskirts of this GC exhibit similar X-ray fluxes. Also, 
using the same comparison as in Sec 13.2.2]. the X-ray to bolometric flux ratio, Fx/F\> Q \ = 9xl0 -5 
is consistent with a main-sequence fore ground star, for which the X-ray to bolometric flux ratio 
saturates at 10" 3 (jvilhu Walterlll987h . 



An analysis similar to the one performed for XMMU 171411—293159 (§ I3.2() excludes the 
possibility that the system is coronally active. The results of the spectral fits with a 2-temperature 
Raymond-smith (RS) plasma are as follows. For kT 2 = 0.93 keV (xl/dof (prob.) = 1.25/14 (0.23)), 
N H>22 = l-15io:!> kT ^ = 0-nlao6 keV i EM i = 2.4^ 2 1 9 xl0 59 cm- 3 and EM 2 < 3.9xl0 55 ' 



In the case of kT 2 = 3.45 keV (xt/dof (prob.) = 1.25/14 (0.23)), N H ,22 = 1.12 



0.1li°^keV, EMi 



2.1+i 8 5 5 xl0 59 



+1.12 
-0.53' 



cm 



cm 
<10~ 3 



3 and EM2 < 5.5xl0 54 cm 3 . Therefore, the distance 
in both cases and lower than the value EM2/EM1 > 1 



independent value EM2/EM1 is 
observed for known coronally active stars (jDempsev et al.l ll993). Thus, the X-ray spectrum of 
this source is inconsistent with that of typical coronally active stars. In addition, the best fit 
galactic absorption iV# values, Nh,22 = 1-15_ 36 or Nh,22 = 1.12_ ' 53 , for the two cases of the 
two-tempeture RS plasma, are significantly larger than the value in the direction of NGC 6304 
{Nh,22 = 0.266). Therefore, assuming an uniform distribution of absorbing material, the best fit 
value of Njj provides further support that XMMU 171421—292917 is unlikely to be a coronally 
active star in the foreground of the GC. 

We therefore conclude that that XMMU 171421—292917 is not a coronally active foreground 
star. Finally, since the X-ray source is spectrally consistent with a qLMXB at a distance of 6kpc, 
we identify the X-ray source as a candidate qLMXB in NGC 6304 with an identified IR companion 
in NGC 6304. 



To estimate the accretion rate, the series of relations in IVerbunt &: van den Heuvell (|1995l ) 
(X-Ray binaries, section 11.3.3) are used. From the radius of the compani on star (R 2 = 9 -7jBq ); 
calculated with the temperature T e g- = 6915 K (which is worked out from iMasana et al.1 (J2006J)) 
and the luminosity L^oi = 7.96xl0 35 erg s" 1 , the companion core mass: M c = 0.27 Mq (Z = 0.02) 
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and M c = 0.31 M Q (Z = 0.0001) are estimated, where y = In (M c /0.25 M ) with ln(i? 2 /-R©) and 
In (L 2 /L Q ) being third order polynomials of y where the polynomial coefficients are theoretical 
fits given in Iverbunt van den Heuvell (|1995l ). From M c and M c = 1.37xl0 -11 (L/L e ) M yr _1 , 
R2/R2 is estimated giving the evolution time-scale of the companion: r ~ 2xl0 7 yr (Z = 0.02) and 
r 2.3xl0 7 yr (Z = 0.0001). Finally, the following relation: 



^1 - _ 2 — ( - - — 

R 2 ~ M 2 \Q Mi 



(3) 



is used to calculate the accretion rate M 2 from the companion on to the neutron star, assuming M 2 = 
O.8M . We found that M 2 = -4.4xl0" 9 Mayr" 1 (Z = 0.02) and M 2 = -3.8xl0~ 9 Mayr" 1 
(Z = 0.0001). With this value of M 2 (for Z = 0.02), the quie scent luminosity in th e DCH model 
is estimated using models that include core neutrino emission ([Yakovlev et al.ll2003l ). The thermal 



X-ray luminosity, -^thermal = 1.1x10 erg s , is roughly c onsistent with ther mal emission from 



other transients with high time-average mass transfer rates (IHeinke et all 120071 ). 



4. Discussion 
4.1. Comparison with ROSAT Observations 



NGC 6304 was observed once previously, using ROSAT/HRI (IRappaport et all 11994 R94 
hereafter), in which four X-ray sources were discovered. Three are spatially coincident with three 
X-ray sources detected in the present observations (in fact, the three correspond to the candidate 
qLMXBs). The fourth X-ray source in R94, however, appears to have faded significantly; the 
present observation detects no X-ray source consistent in position with ROSAT source D. The 
following analysis compares the position and fluxes of the for ROSAT/HRI sources. 

The positions of the four sources in R94 have been examined and compared to those of X-ray 
sources in the present work. No positional uncertainties are given from analysis of R94. Using the 
archived ROSAT/HRI data at HEASARC, the positional uncertainties required for the comparison 
with the present XMM-Newton data were obtained. The source positions are consistent between 
the two different observations: source A, B and C in R94 corresponding to sources #4, #5 and #9 
in this work. The boresight systematic uncertainty (6 arcsec) have been taken into account. For 
source D, no corresponding source is detected within 1 arcmin of the position for source D given by 
R94. A total of 304 counts are found in a 25 arcsec area around the position given by R94, while 
an average of 294.8±2.2 are due to background using a nearby off-source area. This leaves 9±18 
counts due to a possible X-ray source at this position, consistent with no source detection. 

There is a 3cr upper limit on the flux from source D of Fx < 2.4xl0 _14 erg cm -2 s _1 (0.5- 
2.4 keV), assuming a thermal bremsstrahlung spectrum with a temperature kT = 3keV (R94). By 
inspection of the MOS2 data, there is no evidence of any X-ray source in the vicinity of the position 
of source D. 
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Table [5] compares the predicted number of source counts which should be detected with 
ROSAT/HRl during the 5030 sec observation of R94, assuming the observed XMM/pn-med coun- 
trates and the XMM/pn source spectra, including both countrate and spectral uncertainties. For 
sources #4, #5, and #9, the expected number of counts (20.7±1.5, 3.9±0.9, 11.4±1.0, respectively) 
are consistent, within Poisson uncertainties, with the number of source counts observed by R94 (15, 
7, and 17 respectively). 

For source D, which is not detected, the 3a upper limit on the number of counts due to an X-ray 
source at this position (<63 counts, in 11037 sec of integration), implies that <1.8 counts should 
have been detected with ROSAT/HRl, whereas the observed number of counts with ROSAT/HRl 
correspond to 14 count, at the HRI field centre]. In art absence of an error analysis on the 
number of detected source counts from R94, it is not possible to correctly estimate the amount of 
fading; however, taking the detected number of counts at face value implies that source D faded in 
luminosity by a factor of ~10. 

In conclusion, the three qLMXBs are consistent with having the same luminosity between the 
ROSAT/HRl observation and the XMM/pn observation; and source D appears to have faded by a 
factor of ~10. 



4.2. Expected number of quiescent Neutron Stars in Globular Clusters 



The number of X-ray binaries in a GC can be predicted from its internal properties. The 
encounter rate for a GC is pred icted by r oc p\/v where po is the central luminosity density 
and v is the velocity dispersion (jVerbuntl 12002 ) . For a virialized cluster, this relation simplifies 
to r oc Po r c w here r c is the core radius (in physical units, not angular distance). The number 
of quiescent NS s found in NGC 6304 is compared to th e expected number of such objects from 



previous works (jGendre et al. 



2003a 



Heinke et al.ll2003d ). 



In the first one (Table 2 and Figure 3 in iGendre et al.l (j2003al )). the authors found a linear 
relation between the number of quiescent NSs and the encounter rate: N q Ns ~ 0.04 x T + 0.2. Using 
the same normalisation (Fngc 6440 = 100), the encounter rate of NGC 6304 is Tngc 6304 = 6.4. 
The predicted number of qLMXBs is therefore about 0.46. Assuming that three qLMXBs were 
found, this number seems slightly in excess of the prediction but the probability of finding such 
objects when the average expected number is 0.46 is 3.2 per cent , assuming Poisson statistics. This 
probability is not small enough to affirm that there is an inconsistency between the expected and 
the actual number. 



The number of counts given by R94 are corrected for scattering and vignetting by a factor >1, but which is not 
given by R94. Therefore while R94 gives the number of counts which would have been detected from this source if it 
had been located at the centre of the FOV, the number of counts detected was not noted, and so the uncertainty in 
the number of counts has not been determined. 
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The second paper (Table 1 in lHeinke et al.l (|2003d )) compared the encounter rate (normalised 
to the galactic value) to the number of accreting NS systems. The linear relation obtained from 
the data collected in this previous work is simply N = 0.993 x T — 0.046, for which the expected 
number of accreting NSs in NGC 6304 is 0.38. Again, three qLMXBs is a number that seems larger 
than the expectation, but certainly not inconsistent. 

Finally, the en counter rate can be compared to the number of X-ray sources in the GC 
(jPoolev et al.ll2003l ). The expected number of X-ray sources for an encounter rate Tngc 6304 = 45 
is about 9.5, comparable with the seventeen sources detected in the FOV. The following subsection 
discuss the number of background sources that are expected to be detected in the field of view of 
the observation. 



4.3. Background X-ray sources expected in the field of NGC 6304 

To estimate the number of expected Active Galactic Nuclei (AGN) detected during the obser- 
vation of NGC 6304, the field of the pn camera is divided into 3 concentric annulli centred at the 
centre of the exposure map (with radii between and 5, 5 and 10, and between 10 and 13 arcmin). 
The limiting flux in each annulus is estimated using the lowest count rate detected (4.3 counts per 
sec) and the tool webPIMMS, to obtain the absorbed flux a power-law of photon index 2. The 
limiting fluxes are 1.8xl0 _14 erg cm~ 2 s _1 , 1.2xl0 _14 erg cm -2 s _1 , and 0.9xl0 _14 erg cm -2 s _1 , 
for the three annulli respect ively. Using the number-flux (LogN-LogS) distribution function for 
X-ray AGN (jrlasinger 1993 . H93 hereafter), the number of expected background X-ray sources 



is 0.87±0.1, 5.22±0.28 and 10.7±0.75 in the three annulli, respectively, corresponding to a total 
of 17±1 over the whole field. This result is consistent with the number of sources detected in the 
field, and with the tentative conclusion claiming that 3 qLMXBs are detected in the GC NGC 6304. 
However, the error bars only takes into account the uncertainty in the LogN-LogS relation and not 
in the limiting fluxes and are therefore likely to be underestimated. 



5. Summary and Conclusions 

The globular cluster NGC 6304 hosts three candidate qLMXBs. The neutron star radii and 
temperatures obtained for the three candidates are consistent with typical values for neutron stars 
and in accordance with radii of other known field and GCs qLMXBs. No variability is measured 
from any of the candidates over the time scale of the observation, with very weak limits; the large X- 
ray variability in the XMM-Newton /pn background precludes a more detailed variability analysis. 
A comparison with ROSAT observations 14 yr prior (R94), finds the fluxes of the candidates are 
consistent with those of the earlier observation (R94). The absence of a detailed error analysis in 
the fluxes from the published ROSAT analysis precludes placing statistical limit on the amount 
of variability; however, as the number of detected counts in the ROSAT observations was small 
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(<20), the magnitude of intensity variability is limited to be less than a factor of ~few on the 14-yr 
time-scale. 

XMMU 171433-292747 is located in the core (0.79 r c ) of the GC. Its spectrum is acceptably 
fitted with a NS Hydrogen atmosphere model with kT e g = 122^^ eV and i?oo = 11. 6^40 km, 
combined with a power-law component of photon index a = 1.2^Qg. The power-law component 
dominates the spectrum at photon energies above 2keV, and represents 49 per cent of the observed 
flux (0.5-10 keV). This contribution, although its source is not perfectly understood, is similar to 
that found in some other qLM XBs, and is comparable in the fraction of the flux for which it accounts 



to the field qLMXB C en X-4 (jAsai et alJll996bl : lRutledge et alJl2001bl : iMenou k McClintockll2001 



Campana et al.l 12004). While the high stellar density in the core of NGC 6304 precludes identifi- 
cation of a 2MASS counterpart for XMMU 171433-292747, a visual overlap with a 2MASS image 
of the GC shows a possible association with an uncatalogued - and likely unresolved - star. The 
association will require further investigation with higher resolution optical/IR imaging. 

XMMU 171421—292917, a second candidate, has a low signal-to-noise and correspondingly 
larger error bars for the X-ray spectral parameters used to identify it as a qLMXB candidate: 
the NS H-atmosphere temperature is 701^ e V, and i2oo = 23^^ km. A faint 2MASS counterpart 
(mj = 13.063) is identified (with 99.916 per cent confidence), and the X-ray emission is unlikely to 
emerge from the giant star itself. 



XMMU 171411—293159, a third candidate, is a bright X-ray source consistent with a neutron 
star H atmosphere spectral model at the distance of the GC. Its measured properties are A;T e fj = 
U5±HeV and = 10.7±|;f km. The 2MASS IR counterpart (with 99.988 per cent confidence) 
appears to belong to the post-asymptotic giant branch; this is the first qLMXB with this type 
of companion. Other explanations for the source are investigated, concluding that a coronally 
active field star is not a possible explanation, due to the unusual optical/IR colors and unusual 
X-ray spectrum. Optical/IR spectroscopy of 2MASS 17141152—2931594 can definitively determine 
whether the star is a foreground main sequence star, or an evolved post-AGB star at the distance 
of NGC 6304. 



A faint unidentified X-ray source showing an unusual excess of high energy photons resulting 
in a photon index of a = — 2.0^ 2 1- ^ ne nature of this source is unclear. There is no evidence the 
X-ray source was variable during the observation. No off-band counterparts in the 2MASS or DSS 
catalogues are found within 3a of the X-ray position. Deep observations in the off-bands, as well 
as repeated detection and study in the X-ray, can inform interpretation of this source. 

The qLMXB candidates in NGC 6304 are astronomically interesting. First, the number of 
candidates is consistent with predictions from the calculated encounter rate, suppor ting previous 



observational analyses which characterised the qLMXB rate vs. encounter rate in GCs (jPoolev et al 



2003; lHeinke et al.ll2003d : iGendre et al.ll2003al ). This indicates that the qLMXB rate vs. encounter 
rate is not merely a parameterisation resulting from observations, but can predict the number of 
qLMXBs in unobserved GCs. 
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NGC 6304 contains one of the most luminous known qLMXBs. XMMU 171433-292747, with a 
intrinsic luminosity Lx = 10 33 erg s _1 (0.5-10 keV), is the fourth most luminous GC qLMXB after 
47 TucX7 (L x = 1.5xl0 33 erg s^ 1 (0.5-10 keV)), 47 Tuc X5 (L x = 1.4xl0 33 erg s" 1 (0.5-2.5 keV)) 
and the qLMXB in M28 (L x = 1.2xl0 33 erg s" 1 (0.5-10 keV)). These represent the upper limit 
of the expected luminosity for LMXB in quiescence. In addition, XMMU 171433—292747, with its 
high flux (F x = 2.3 xl0~ 13 erg cnT 2 s~\ after 47 Tuc X7, 47 Tuc X5 and the qLMXB in M28 with 
fluxes 5.3, 4.3 and 3.3^{'^ xl0~ 13 erg cm -2 s -1 , respectively - see Tabled]), is a good candidate to 
constrain the equation of state of dense matter. 

NGC 6304 hosts two candidate qLMXBs lying at large distance from the hosting core, but still 
within the tidal radius; XMMU 171411-293159 and XMMU 171421-292917 are located at ~ 32r c 
and ~ 15r c , respectively. Previously, all GC qLMXBs were located within < 7r c of the optical 
centre of their GC (see Table [6]), the most distant belonging to U24 in NGC 6397, located 6.8r c 
from the centre; noting, however, that NGC 6397 is a core collapse GC with r c = 0.05'. To sum 
up, the qLMXB candidates in NGC 6304 are of much greater distance from the GC centre than 
has been observed previously, and may require revisiting theory of qLMXB formation and binary 
evolution in GC. 

Deeper X-ray exposures are required to confirm the classification of the candidates qLMXB 
found in NGC 6304. Specifically, it is necessary to constrain to an uncertainty of few percent 
(as already performed for confirmed qLMXBs in to Cen and 47 Tuc, see Table [4]) rather than few 
tens of percent here. For XMMU 171433—292747, better angular resolution will confirm the iso- 
lated nature of the candidate qLMXB, or may resolve multiple X-ray sources in the core. New deep 
observations of XMMU 171421-292917 and XMMU 171411-293159 will produce a localisation 
unbiased by bad detector columns. At optical wavelengths, observations of the aforementioned 
candidates can determine unequivocally the evolutionary state of the counterparts (a giant for 
2MASS 17142095-2929163, and a post-AGB star for 2MASS 17141152-2931594), as well as mea- 
sure their binary orbital periods through ellipsoidal variations. The spatial isolation and brightness 
of the IR counterparts implies this work can be done easily using ground based, moderate sized 
telescopes. Optical spectroscopic observations can confirm that the counterparts have radial ve- 
locities consistent with that of NGC 6304. Finally, high spatial resolution observations of the GC 
core in the optical/IR bands can pinpoint a possible counterpart for XMMU 171433—292747 and 
confirm the qLMXB nature of this X-ray source. 
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Fig. 1. — XMM/pn image of the globular cluster NGC 6304. Each small circle shows the position of 
the 17 detected sources. The small circles labelled Ml through M6 correspond to sources detected 
on the MOS cameras only. The bigger solid circle and large dashed circle represent the half mass 
radius and the tidal radius (respectively) of NGC 6304. The v alues thm = 1-41' and r t = 13.25' 
are from the Harris catalogue of globular clusters ([Harris! Il996l ). The source numbers correspond 
to those in Table Q] and Table [2j 
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Fig. 2.— Folded spectrum of source #8 (XMMU 171413-293415) fitted with an absorbed power 
law model (a = — 2.0^22) with the H-absorption fixed Nh,22 = 0.266. The binning of the counts 
is different from the standard binning (minimum of 20 or 25 counts per bin). Due to the low 
signal-to-noise of this source, bins of 1 keV width (except for the 0.5 to 1 keV bin) are chosen. 
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Fig. 3. — Lightcurve between 5 and lOkeV, with 100 sec bins. Bottom: lightcurve of extraction cir- 
cular region of source #8 (25 arcsec). The horizontal line at 0.0121 counts/sec (0.0078 background 
counts/sec expected) is the mean count rate. Top: lightcurve of the background area chosen for 
source #8. 
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Fig. 4.— Folded model to the pn camera spectrum of source #4 (XMMU 171433-292747). The 
solid line is the best fit model of a neutron star atmosphere nsa (Nh,22 = 0.266 held fixed) together 
with a power-law component of slope a = 1.2+Qg. Due to the low signal-to-noise at large energies 
(above 2.5 keV), the last 9 bins are grouped into three separate bins. 
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Fig. 5. — Folded model to the pn camera spectrum of source #9 (XMMU 171411—293159). The 
solid line is the best fit (tabulated) model of a NS H-atmosphere (Njj,22 = 0.266 held fixed). Due 
to the low signal-to-noise at large energies, the last three bins are grouped into a single bin above 
4keV. 
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Fig. 6. — Color (J-K) Magnitude (K-band) diagram of the globular cluster NGC 6304 showing 



isoch rones (small squares) for ages in the range t = 10 8,55 — 10 1015 yr and Z = 0.00488 (jMarigo et al 



2008). All 2MASS catalogued stars (small dots) within 6.5 arcmin of the optical centre of the 



cluster are included, down to a magnitude of mx ~ 15. The large squares with error bars are 
the IR counterparts to source #9 and #5, with mj = 8.769 and mx = 13.063, respectively. This 
diagram suggests that the counterpart for XMMU 171411—293159 belongs to the post-asymptotic 
giant branch. A spatial analysis of the stars that lie redward of the giant branch compared to those 
that lie blueward shows that there is an over-density of bluer stars near the GC core. In other 
words, redder stars are less likely to belong to the GC and are likely to be contaminating field 
stars. Thus, contamination by field stars can explain the poor matching between the theoretical 
isochrones and the actual data. 
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Fig. 7. — Color (B-V) Magnitude (V-band) diagram of the globular cluster NGC 6304 with the 
isoch rones (small squares) for ages in the range t = io 8 - 55 -10 10 - 15 yr and Z = 0.00488 dMarigo et al 



2008). Also included are all the stars (s mall dots) within 6.5 arcmin of the optical centre (S. 



Ortolani, private communication, 2008 and lOrtolani et al.1 (|2000j)) down to a limiting magnitude of 

V ~ 21. The large squares with error bars are the optical counterparts to source #9 and #5, with 

V = 10.65 and V = 14.85, respectively. Similarly to FigEl this suggests that the counterpart for 
XMMU 171411-293159 is a post-AGB star at the distance of NGC 6304. Also, field contamination 
could be responsible for the poor matching between the theoretical isochrones and the actual 2MASS 
data. 
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Fig. 8. — Color (J — K) Color (B — V) diagram showing isochrones of low-mass (M < 5M Q ) 
stars at t = 10 7 yr (small points) for a range of metallicities (Z = 0.0001 — 0.03, A.Z = 
0.0005). The large points with error bars are corresponding to XMMU 171421—292917 (#5) and 
XMMU 171411—293159 (#9). In considering source #9 as a foreground coronally active star, it lies 
marginally off the main sequence (between 2-3<r) at zero reddening (E(B — V) =0). The discrep- 
ancy worsens if source #9 is at greater reddening; the vector beginning at source #9 terminates at 
the intrinsic color 2MASS 17141152-2931594 would be at E(B - V) = 0.53, that of NGC 6304. 
This does not support the hypothesis that XMMU 171411—293159 is a typical coronally active 
star, in the foreground of NGC 6304. 
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Fig. 9.— Folded model Spectrum of source #5 (XMMU 171421-292917). The solid line is the 
best fit (tabulated) model of a neutron star atmosphere {Nh, 22 = 0.266 held fixed). 
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Table 1. X-ray Catalogue of sources in the field of NGC 6304 (Astrometry) 



Object Name 


R.A. 
(J2000) 


Decl. 
(J2000) 


<5R,.A.\5Docl. a 
(arcsec) 


A/r c b 


ID 


XMMU 171403- 


-292438 


258.51323 


-29.41030 


±0.5\0.5 


33.4 


10 


XMMU 171411- 


-293159 


258.54803 


-29.53317 


±0.5\0.4 


31.8 


9 


XMMU 171413- 


-293415 


258.55432 


-29.57151 


±1.0\0.8 


38.6 


8 


XMMU 171415- 


-292232 


258.56320 


-29.37568 


±0.9\0.8 


31.9 


7 


XMMU 171417- 


-293222 


258.57100 


-29.53968 


±0.9\1.0 


28.5 


6 


XMMU 171421- 


-292917 


258.58731 


-29.48801 


±0.8\1.0 


15.2 


5 


XMMU 171433- 


-292747 


258.63636 


-29.46310 


±0.7\0.6 


0.79 


4 


XMMU 171449- 


-292001 


258.70731 


-29.33378 


±1.0\0.9 


42.3 


3 


XMMU 171451- 


-293653 


258.71296 


-29.61490 


±1.0\0.9 


49.1 


2 


XMMU 171453- 


-292436 


258.72141 


-29.41001 


±0.7\0.6 


29.2 


1 


XMMU 171516- 


-292224 


258.81800 


-29.37320 


±0.5\0.5 


52.2 






Note. — The sources are sorted by increasing R.A. Both coordinates (R.A. and 
Decl.) are given in decimal degrees, after the astrometric correction. 

Statistical uncertainty from the source detection in seconds of arc. 

b distance from GC centre in r c units - r c =0.21' (Harris 96). 



Table 2. X-ray Catalogue of sources in the field of NGC 6304 detected with the MOS camera 

only 



Object Name 


R.A. 


Decl. 


<5R.A.\ 5 Docl. a 


S/N 


Rate 


ID 






(J2000) 


(J2000) 


(arcsec) 




(cts ks- 1 ) 




XMMU 171514- 


-293254 


258.81098 


-29.54841 


±1.1\0.8 


4.3 


2.8±0.8 


Ml 


XMMU 171418- 


-293053 


258.57689 


-29.51474 


±0.7\0.8 


7.1 


3.1±0.6 


M2 


XMMU 171443- 


-292511 


258.68223 


-29.41987 


±0.6\0.6 


5.1 


2.0±0.5 


M3 


XMMU 171401- 


-292145 


258.50750 


-29.36265 


±1.0\0.6 


6.3 


3.3±0.7 


M4 


XMMU 171356- 


-291720 


258.48705 


-29.28889 


±0.7\1.0 


4.5 


3.4±0.9 


M5 


XMMU 171505- 


-292235 


258.77451 


-29.37666 


±1.2\0.7 


5.6 


3.6±0.8 


M6 



Note. — The right ascension (R.A.) and declination (Decl.) are given in decimal degrees. 
Statistical uncertainty from the source detection in seconds of arc. 
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Table 3. X-ray Catalogue of sources in the field of NGC 6304 (Fluxes) 



ID 


S/N 


Rate a 
(ctsks- 1 ) 


kT cB b 
(eV) 




xS/dof (prob.) 


Fx c 


Type d 


10 


5.2 


4.5±1.9 




(1) * 


0.74/16 (0.76) 


0.37 


?? 


9 


20.2 


25.2±2.2 






1.20/31 (0.21) 


1.52 


qNS? 


8 


5.3 


7.1±2.0 




-a.o±|:i 


1.73/7 (0.10) 


2.8 


?? 


7 


4.3 


4.3±1.9 




(i) * 


0.84/16 (0.63) 


0.37 


?? 


6 


8.0 


7.2±1.7 




(i) * 


0.75/13 (0.71) 


0.76 


?? 


5 


7.0 


9.1±2.0 


70t» 




1.09/16 (0.36) 


0.59 


qNS? 


4 


19.3 


43.8±2.8 


1 99+31 


1 2+°' 7 


0.85/42 (0.75) 


2.3 


qNS? 


3 


4.3 


6.2±1.9 




(1) * 


0.68/16 (0.82) 


0.77 


?? 


2 


5.4 


6.1±1.9 




(1) * 


1.37/15 (0.15) 


0.86 


?? 


1 


8.3 


15.3±2.1 




0.52±0.34 


0.96/15 (0.50) 


2.47 


?? 





4.7 


10.3±1.7 




-0 3+ - 7 


0.84/12 (0.60) 


4.2 


?? 



Note. — Nh column density fixed to Nu,22 = 0.266 for all spectral fits. 
a Background subtracted count rate in counts per ks. 

b Bcst fit parameters. kT e g for a neutron star atmosphere model and a for a power-law 
component. Both were combined in the case of source #4. 

c Unabsorbcd flux (units of 10~ 13 erg cm~ 2 s _1 ) between 0.5-10 kcV. 

d Type of X-ray source. qNS = transient neutron star in quiescence. 

*The value of the slope in parenthesis arc held fixed at the value a = 1. 



-38- 



Table 4. H-atmosphere spectral parameters for source #4, #9 and #5 for other known probable 

GC and field qNSs 



Object 


Roo 

(km) 


(eV) 


d 

(kpc) 


Nh, 22 


a 


Lx 


Fx 




Rcf. 


YMA/TT7 171411 9Q31 
-a.iv.ij.vi \j l i l*±± ± — zyoioy 


iU -'-3.1 


iic+21 
ii0 -16 


^ Q7 

o.y / 


(a 




u.uo 


1 52 


l 


present work 


YMMTT 171491 9Q9Q1 7 
yvivuvi \j 1 1 i.'-tZi i — zyzy ± i 


/ci -10 


7n +28 


^ Q7 






n ri7 

U.U I 


u.oy 


l 


present work 


YA/TA/TTT 1 71 /1QQ 0007/17 

A1V11V1U 1 ( 14oo — zyz (4 / 


n - b -4.6 


U2 -45 


c; ti7 
o.y / 


(O.ZDD j 


i 9+0.7 
L2 -0.8 


1.U1 


z.o 


l 


present work 




14.3_ 2 x 


fifi+ 4 
66_ 5 





^u.uy J 




n tr _i_n o 


1.0 ( itU.D / 


2 


i 
i 


r,jCVn('F<PTC i 't 


13 6+°- 3 


g 7 +2 


5.3 






n ^94-0 09 


u . yoiu.uu 


2 


2 


IVilo ^-CjJr 1^ J 


12.8_ 4 


7fi +3 
3 


7 7 


/n m 1^ 




U. < oztU.UD 


l.UoztU.Uo 


2 


q 
o 


47 Tuc X5 


19.0 + 7g 




4.85 


0.09+°°5 




1.4 


4.3 


1,* 


4 


47 Tuc X7 




105.4l 5 5 .6 


4.85 


042+ ' 018 




1.5 


5.3 


1:* 


5 


M 28 (#26) 


14.5«J 


gn+30 


5.5 


0.26+.0.04 




1 2+°' 7 


3.3511;? 


3 


6 


M 30 A-l 


13.41" 


94+17 


9.0±0.5 


029+ - 017 
u.uza_ Q 012 




0.71 


0.73 


1:* 


7 


NGC 6397 (U24) 


4.9+i 4 


57-92 


2.5 


0.1-0.26 




0.08 


1.06 


5,* 


8 


M80 CX2 


(10) 


82±2 


10.3+°-? 


0.09l° 025 




0.29±0.02 


0.23±0.02 


4 


9 


M80 CX6 


(10) 


76l 5 6 


10.31°!? 


n 99+O.O8 
u -^-0.()7 




0.09±0.01 


0.07±0.01 


4 


9 


NGC 2808 C2 


(12) 


81.3l|| 


9.6 


0.82±0.40 




0.26±0.04 


0.24±0.04 


1 


10 


NGC 3201 16 


(10) 


■i 7o + - ;) ° 
i ' Z -42 


5.0 


(1.17) 


1.04±0.63 


0.34 


1.0 


0,* 


11 


Aql X-l 


13.4+° 


135li« 


5 


u.oo_ 07 


(1.0) 


4.4 


14.7 


1.* 


12 


Cen X-4 


16.8l^ 


59l 6 6 


1.2 


(0.055) 


1 0+°' 6 


0.17 


0.6 


1,* 


13 


a 4U 2129+47 


y . d -4.6 


60+iJ 


1.5 


(0.28) 




o.07l°;J 5 3 


2 4+5-1 
z -^-1.0 


6 


14 


a 4U 2129+47 


16±^ 4 g 


soil 9 . 


6.0 


(0.17) 




0.56i°i 2 


i n+0.3 
1 — i -0.T 


6 


14 


a 4U 1608-522 


12-31I1 


1301^° 


3.6 


(0.8) 




0.83±0.42 


5.4±2.7 


1 


15 


a 4U 1608-522 


651SJ 


solil 


3.6 


(1.5) 




0.73±0.47 


4.7±3.0 


1 


15 



Note. — For the qLMXBs of the present work, the mass was kept fixed to 1.4 Mq. As mentioned in the text, the tabulated model of 
the neutron star atmosphere was used for XMMU 171411-293159 and XMMU 171421-292917 while a nsa+PL model was used to fit the 
spectrum of XMMU 171433—292747. The top part of the table contains GC qLMXBs while the second part only shows the parameters 
of field qLMXBs. The inferred radius is measured using the distance to the GC quoted in the text. Previously measured values of the 
radius have been converted from R to Roo, the projected radius. Error bars are 90 per cent confidence. Values in parenthesis are held 
fixed. The unabsorbed X-ray luminosity Lx and flux Fx are expressed in units of 10 33 erg s _1 and 10 — 13 erg cm~ 2 s — 1 , respectively, 
both in X-ray bands listed below. 

a Diffcrcnt assumed values of D and/or Njj. 

°0.2-10kcV band. 

^.S-lOkcV band. 

2 0.1-5kcV band. 

3 0.5-8keV band. 

4 0.5-6keV band. 

5 0.5-2.5keV band. 

6 0.5-3.0keV band. 

*No error was provided by the authors. 
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v Indicatcs that the source flux varies. 



Note. — Re f erence s: 1 iRutledee et al.l [|2002al); 2, iGendre et all |2003t ); 3. iGendre et alj J2()03ah ; 4. iHeinke et alj d2003bh 
5, IHeinke et al.l 6. jBecker_etal1 j2003h ; 7 . iLuegeT et al 



, .Grindlav et alj ll2001bh; 9. IHeinke et alj j2003 j); ~10 , 
13. iRutledge et al.l j2001bll ; 14. IRutledee et alj J200cl'> ; 



(2007 



ISeryillat et alj iboOSt); 1 1, IWebb et alj I2OO6I) 12. IRutledge et~aU j2001 
15. IRutledee et alj jl99sh 
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Table 5. Difference in the Luminosities between the RO SAT and the current XMM-Newton 

observations 



XMM-pn ROSAT 



ID 


fcT cffiBB a 


ECF b 


IDr94 


Counts predicted 


Counts Observed d 


4 


260±30 


0.0935±0.0015 


A 


20.7±1.5 


15 


5 


170±30 


0.0855±0.0065 


B 


3.9±0.9 


7 


9 


200±30 


0.090±0.002 


C 


11.4±1.0 


17 






0.062 


D 


<1.8 


14 



Note. — For source D, due to the non detection, a 3<r upper limit is calculated. 

a The fcT e ff BB is the black-body effective temperature, from spectral fits to the present 
XMM data. 

b The energy-conversion-factor (ECF) is the number of ROSAT/HRI counts per XMM 
count for each source, which depends on fcT c ff BB an d Nh (assumed to be Nh,22 = 
0.266). For R94 source D, undetected in the present observation, an absorbed thermal 
bremsstrahlung spectrum with kT = 3 keV is assumed. 

c The total counts predicted for the R94 observation based on the present XMM ob- 
servation, using the source spectrum measured here; the predicted HRI countrate from 
WebPIMMS; and the ROSAT/HRI observation duration of 5030 sec. 

d The total number of source counts from ROSAT/HRI observation (R94), corrected for 
observation at the HRI centre (that is, after subtracting background counts, correcting for 
vignetting, quantum efficiency and scattering). The uncertainty in this number is not given 
by R94. 
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Table 6. Distance from the core of the candidate qLMXBs and other known or probable qNSs in 

GCs. 



Object Name 


d c r c A/r c Ref. 




(arcmin) (arcmin) 



XMMU 171411-293159 


6.7 


0.21 


31.8 





XMMU 171433-292747 


0.17 


0.21 


0.79 





XMMU 171421-292917 


3.2 


0.21 


15.2 





wCen 


4.38 


2.58 


1.7 


1 


M13 


0.74 


0.78 


0.95 


2,10 


47 Tuc X7 


< 0.37 


0.38 


< 1.5 


3 


47 Tuc X5 


< 0.37 


0.38 


< 1.5 


3 


M 28 (#26) 


0.05 


0.24 


0.21 


4,10 


M 30 A-l 


0.03 


0.06 


0.5 


5,10 


NGC 6397 (U24) cc 


0.34 


0.05 


6.8 


6,10 


M80 CX2 


0.063 


0.15 


0.42 


7,10 


M80 CX6 


0.324 


0.15 


2.2 


7,10 


NGC 2808 C2 


0.08 


0.26 


0.31 


8,10 


NGC 3201 16 


5.8 


1.43 


4.1 


9,10 



Note. — This table is a compilation of the positional information of 
known qLMXBs and candidates qLMXBs in GCs. The columns, from 
left to right, are the object name (or the GC hosting the qLMXB), 
the distance of the qLMXB from the optical centre of the GC (in 
arcminutes), the core radius of the GC (in arcminutes) and the distance 
of the qLMXB from the optical centre of the GC in units of core radii. 
References : 0, pres ent work; 1 , iRutledge et al.1 d2002all; 2 , 
. IHeinke et al.l j2003b|l; 4,lBecker et al.N2003h; 
Grindlav et al.l ibOOlbh; 7, iHeinke et alJ 
9, IWebb et afl(l2006l) ; 10. lHarrisI 



iGendre et all [12003 



Lugger et al.1 [|2007r); 6 



1 2003311 ; 8. IServillat et"afl fcOOSft ; 
1 1996), for the core radii values. 



C NGC 6397 is a core collapsed star cluster. 



